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ABSTRACT 

We investigate the metallicity properties of host galaxies of long Gamma-ray Bursts 
(GRBs) in the light of the Fundamental Metallicity Relation (FMR), the tight de- 
pendence of metallicity on mass and SFR recently discovered for SDSS galaxies with 
stellar masses above 10 9 ' 2 M Q . As most of the GRB hosts have masses below this 
limit, the FMR can only be used after an extension towards lower masses. At this 
aim, we study the FMR for galaxies with masses down to ~ 1O 8 3 M0, finding that the 
FMR does extend smoothly at lower masses, albeit with a much larger scatter. We 
then compare the resulting FMR with the metallicity properties of 18 host galaxies 
of long GRBs. While the GRB hosts show a systematic offset with respect to the 
mass-metallicity relation, they are fully consistent with the FMR. This shows that the 
difference with the mass-metallicity relation is due to higher than average SFRs, and 
that GRBs with optical afterglows do not preferentially select low-metallicity hosts 
among the star-forming galaxies. The apparent low metallicity is therefore a conse- 
quence of the occurrence of long GRB in low mass, actively star-forming galaxies, 
known to dominate the current cosmic SFR. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic explo- 
sions in the Universe (see Zhang & Meszaros 2004 for a re- 
view) and detected in the 7-rays with a frequency of about 
one per day over the whole sky. The 7— ray emission is ac- 
companied by a long-lasting tail, called afterglow, usually 
detected over the whole electromagnetic spectrum. Their 
extreme brightness easily over-shine the luminosity of their 
host galaxy and makes them detectable up to extreme high 
redshift as shown by the discovery of GRB 090423 at z = 8.2 
( Salvaterra et al.|2009| |Tanvir et al.|2009[). GRBs are u sually 



divided in two, broad classes ( Kouveliotou et al.|1993[ ): short 
GRBs, which are believed to result from the merger of two 
compact objects, and long GRBs, associated to the collapse 
of the core of a massive star, as a Wolf Rayet s tar (|Yoon| 



et al. 2006 Woosley fc Hcgcr 2006 Yoon et al. 20081. In 



this paper, we limit our analysis to the class of long GRBs. 

Recent studies on the final evolutionary stages of mas- 
sive stars ( |Woosley fc Heger||2006| |Fryer et aL]|1999[ ) have 
suggested that a Wolf-Rayet star can produce a long GRB if 
its mass loss rate is small, which is possible only if the metal- 
licity of the star is lower than ~ 0.1 — 0.3 Z@. In this case, 



the specific angular momentum of the progenitor allows the 
loss of the hydrogen envelope while preserving the helium 
core. In this view, GRBs may occur preferentially in galax- 
ies with low-metallicity ( Fynbo et al.|2003 Prochaska et al 



2004| |Fruchter et al.|2006| |Stanek et al.|2d06[ >, although we 
have to stress that low-metallicity progenitors do not neces- 
sarily imply low-metallicity host galaxies. Indeed, owing to 
the existence of metallicity gradients inside galaxies, GRBs 
could form from low-metallicity progenitors also in hosts 
with relatively high metallicities ( |Campisi et al.|[2009[ ). 

Up to now, we have been able to detect the host galaxy 
of ~ 70 long GRBs with known redshift. In more than half of 
the cases, the observations allowed to determine the stellar 
mass and the star formation rate of the galaxjQ The ob- 
servational information gathered so far indicates that long 
GRBs with optical afterglows are typically found to re- 
side in low-mass, dwarf galaxies with average stellar masses 
M t ~ 1 — 5 x 10 9 Mq and high specific star formation rate 
(SSFR=SFR/A/*). Information about the chemical content 
of this objects are known only for a subsample of the hosts 
( Savaglio et"ah1|2009| [Levesque et al.||2010d|c[ ). While most 
of the long GRBs are in low-metallicity galaxies, a few cases 
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for which the galaxy metallicity is found to be quite high 



do exist (e.g. GRB 020819, Levesque et al. 2010c 



Yoldas, et al.pOlO and GRB 050401, Watson et al 



Kiipcii 



2006), 



so that the role of metallicity in driving the GRB phenom- 



ena remains unclear and it is still debated (Fynbo et al 



2003||Prochaska et al.|2004||Stanek et al.|2006||Fynbo et al. 
2006 Wolf fc Podsiadlowski|2007| JPrice et al.|2007||Modjaz 
et al'||2008| |Kocevski et al.|2009| |Savaglio et al.||2009[ |Gr~ 
ham et al.|2009b|a||Levesque et al.|2010a|c|e||Svensson et al.| 
2010||Fan et al.||2010| ). " 



Many recent studies have attempted to find similarities 
and differences between the GRB host population and the 
normal field galaxy one (see, for example, |Fynbo et al.| 
20081. In particular, these studies compared the observed 
mass-metallicity relation (or luminosity-metallicity re- 
lation) of the two populations obtaining contradictory 
results. From the analysis of a whole sample of known GRB 



hosts, Savaglio et al. ( 2009 \ concluded that there is no 



clear indication that GRB host galaxies belong to a special 
population. Their properties are those expected for normal 
star-forming galaxies, from the local to the most distant 
universe. On the other hand, the study of sub-samples with 



well-determined chemical properties (e.g. Levesque et al 
2010a|c Han et al. 2010 1 suggests that most of the long 
GRB host galaxies fall below the M — Z relation for the 
normal galaxy population. 



The aim of this work is to further test the differences 
between GRB hosts and field galaxies by taking advantage 
of the new Fundamental Metallicity Relation (FMR) re- 
cently introduced by |Mannucci et al.|p010[ ). The FMR is a 
tight relation between stellar mass M*, SFR and gas-phase 
metallicity. Local SDSS galaxies define a surface in the 3D 
space of these three quantities, with metallicity well deter- 
mined by stellar mass and SFR. The residual metallicity 
scatter around this surface is very small, about 0.05 dex, 
similar to the expected uncertainties. Also, the same FMR 
defined locally by SDSS galaxies is found to describe, with- 
out any evolution, the properties of high-redshift galaxies, 
up to z=2.5. The origin of the strong, monotonic evolution 
of the mass-metallicity relation over the same redshift range 
(e.g., |Tremonti et al.|[2004[|Erb et al.||2006b| > is due to the 
increase of target SFR with redshift, resulting in sampling 
different parts of the same FMR at different redshifts. At 
even higher redshifts, galaxies are found to evolve off the 
FMR ( jMaiolino et al.|2008| jMannucci et al.|2009[ ), and this 
effect is under test with a larger number of observations 
(Troncoso et al., in prep.). 



The ranges of mass and SFR over which the FMR was 
measured are limited by the number of galaxies in the SDSS- 
DR7 sample used, which become rare at log(M*/M0) be- 
low 9.2 and above 11.4, and at log(SFR/M yr" 1 ) below 
-1.4 and above +0.8. For this reason, a comparison with the 
hosts of GRBs can only be done by extending this relation 
using lower mass galaxies, while a simple extrapolation of 
the FMR of massive galaxies could produce spurious effects. 



EXTENDING THE FMR TOWARDS 
SMALLER MASSES 



To derive the FMR, |Mannucci et al.| ( |2010[ ) have split 
~140000 SDSS-DR7 galaxies into bins of mass and SFR 
having width of 0.15dex in both quantities. To have a good 
estimate of both median and dispersion of the metallicity 
for each value of mass and SFR, only bins containing more 
than 50 galaxies have been used. This severely limits the 
range of mass and SFR over which the FMR has been 
measured, even if a significant number of galaxies outside 
these ranges are present in the original sample. Among 
the ~ 140000 SDSS-DR7 galaxies selected by |Mannucci| 
(20101 requiring 0.07<z<0.30 and signal-to-noise 



et al. 



ratio SNR(Hq)>25, about 2000 (1.4%) have masses below 
10 9 2 M Q . Here we intend to use these galaxies to extend 
the measured FMR. 



Mannucci et al. (20101 have introduced the new quan- 



tity [i, a defined as a linear combination of stellar mass and 
SFR: 



[i a = log(M,) — a log(SFR) 



(1) 



and have demonstrated that, for «=0.32, all galaxies at 
z<2.5 show the same dependence of metallicity on /jo.32 and 
the same range of values of /io.32- In other words, the in- 
troduction of this quantity roughly defines a projection of 
the FMR that minimizes the scatter, i.e, corresponds to ob- 
serving the FMR "edge on" . From a physical point of view, 
metallicity is found to increase with mass and decrease with 
SFR, therefore a combination of these two quantities, with 
a negative factor for SFR, is expected, and found, to show a 
better correlation with metallicity. It is worth noticing that 
the dependence of metallicity only on /10.32 is not exact, as 



no part of the FMR is exact ly a plane (sec Fig. 2 of Mannucci 
et al.|20l0 |, nevertheless this is a convenient approximation 



To avoid binning the limited number of galaxies with 
low mass into a large number of classes of mass and SFR, 
we extend the FMR directly by considering the combina- 
tion /iq.32- We consider the ~1400 galaxies in Ma nnucci] 



et al. (20101 sample with 8.3<^o.32<9.4. This is a small 
sample, side by side to the large sample of ~140.000 galaxies 
with larger values of H0.32, and problems with contamination 
are possible. Indeed, while ~1300 galaxies with low ^0.32 
show low metallicities, 86 of them, corresponding to 0.0007 
of the full sample, have large values of metallicities, above 
12+log(0/H)=8.9, with the same distribution of metallic- 
ity of the large population of massive, quiescent, metal-rich 
galaxies. Given the intrinsic uncertainties on mass and SFR, 
these 86 galaxies are likely to be metal-rich galaxies whose 
H0.32 is incorrectly measured and scattered towards low val- 
ues. We remove these galaxies from the sample, and divide 
the remaining ~1300 galaxies in bins of /io.32, and for each 
bin we compute median and standard deviation of metallic- 
ity. For comparison, we also compute the mass-metallicity 
relation considering the ~1700 galaxies with masses between 
10 8 - 3 and 10 9 - 4 M Q . 

The results are shown in figure [T] where low-mass 
galaxies are compared to galaxies of larger M*. The left 
panel shows the mass-metallicity rel ation. At masses abov e 
~ 10 10 M Q this is fully consistent with|Tremonti et al.||2004|, 



while it shows lower values of metallicity and a steeper 
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Figure 1. Left: Mass-Metallicity:, metallicity of low- mass SDSS galaxies (blue open dots with la dispersions) as a function of stellar 
mass. The coloured lines are local SDSS galaxies from Mannucci et al. (2010), color-coded from red to blue according to increasing SFR. 
The black thick line shows the polynomial fit to the mass-metallicity relation in Mannucci et al. (2010). The black dashed line is the 
mass-metallicity relation in Jabran Zahid et al. (2010), transformed to the same metallicity scale. The host galaxies of long GRBs are 
over-plotted (red solid dots, labelled with the GRB date). The red thick line is a linear fit to these GRB host data, with ±lcr bands 
shown in light red. It is clear that GRB host follow a different relation and show systematically lower metallicities. The lower panel 
shows the difference between the metallicity of the GRB hosts and the mass-metallicity relation of SDSS galaxies, showing the systematic 
offset toward lower metallicities. Right: Fundamental metallicity relation: metallicity as a function of £to.32=log(Af*)-0.321og(SFR) in 
solar units. The black solid line is the linear fit to the low-mass SDSS galaxies. For comparison, the black dotted line is the extrapolation 
of the 2nd degree fit to the FMR of the SDSS galaxies as defined in Mannucci et al. (2010) and plotted for SFR=0. The linear fit to the 
GRB host data (red line with the ±lcr band) shows that GRB hosts are fully compatible with the FMR defined by local SDSS galaxies. 
This is also shown in the lower panel, where the metallicity difference of GRB hosts with the FMR is plotted. 



slope at lower masses. This is probably due to the differ- 
ent selections: our requirement of a high signal-to-noise ratio 
(SNR>25, see |Mannucci et al.|2010[ ) in the Ha line preferen- 
tially select galaxies with high SFR and, as a consequence, 
lower metallicity, especially at low stellar masses. 

On the right of figure [I] the FMR is shown. Low mass 
galaxies appear to extend smoothly the FMR, with a linear 
relation between metallicity and /j.o.32- The resulting FMR 
can be described by: 



\2 + log(0/H) 



= 8.90 + 0.37m -0.14s - 

+0.12ms - 0.054s 2 

= 8.93 + 0.51(^0.32-10) 



0.19m 2 

for ^0.32 ^ 9.5 
for /io.32 < 9.5 



where m=log(M*)-10 and s=log(SFR) in solar units. 

It is evident that the intrinsic scatter around the FMR 
increases towards lower values of ^to.32- The residual scatter 



is larger than the expected errors on metallicity, mass and 
SFR, even if the uncertainties on stellar masses from SED 
fitting could increase towards low masses. This increasing 
scatter towards dwarf galaxies is a well known effect proba- 
bly related to a large spread of histories and current levels of 
star formation (|Hunter fc Hoffmanl|1999| |Hunt et al71|2005 



Zhao et al.|2010 l 



3 THE METALLICITY OF THE HOSTS OF 

(2) GRB 

The formation of long GRBs is thought to be related to the 



collapse of a very massive, low metallicity star (Woosley & 
|Heger|2006||Fryer et al.|1999[ ). Thus, it has been argued that 
the occurrence of a long GRB may be linked to a overall low 
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metal content of its host galaxy, making GRB hosts a biased 
galaxy sample with respect to the normal field population. 
In order to check whether this bias exists, we consider the 
properties of the GRB hosts in the light of the observed 
FMR for normal field galaxies. To this extent, we collect all 
the GRB host galaxies at z<l with available observations 
to measure, at the same time, stellar mass, SFR, and gas 
phase metallicity. 

Line fluxes of long GRB hosts have been published 
by several authors (|Savaglio et aT1|2009| |Han et al.||2010| 
Levesque et al.|2010b |. We have used these compilations to 



Table 1. Properties of the GRB hosts 



measure gas-phase metailicities, using the method described 
in |Maiolino et al.| ( |2008[ ) and |Cresci et aT] (|2010[), and ex 



pressing it in the same scale as in Nagao et al 



ley fc Ellison (2008), and Mannucci et al. (20101, where 



(20061, Kew 



solar metallicity is 12+log(0/H)=8.69. Many metallicity 
indicators have been proposed that are based on line ra 
tios (e.g., |Pettini & Pagel||2004| |Nagao et al.||2006| |Kew 



|ley fc Ellison||2008[ ), buT none of them is without prob- 
lem. For example, R23 has two branches, with two dif- 
ferent metailicities associated to each value of R23. Both 
Hq/[NII]A6584 and [Oil] A3727/[NeIII] A3869 have mono- 
tonic variations with metallicity and no dependence on ex- 
tinction, but they include fainter lines, especially for very 
high or very-low metailicities. [OIII]A4958,5007/[OII]A3727 
and [OIII]A4958,5007/[NII]A6584 are sensitive to extinc- 
tion, which is usually poorly known. Following |Nagao et al.| 
(2006), we measure metailicities by simultaneously consid- 
ering all the flux ratios among the relevant emission lines, 
fitting these values with two free parameters, metallicity and 
extinction. Usually this method can obtain a reliable value 
of metallicity, avoiding or mitigating the intrinsic problems 
of each individual line ratio. In contrast, extinction is usu- 
ally very poorly constrained, because most of the line ratios 
used involve line with similar wavelengths. For this reason, 
when flux ratios between different hydrogen Balmer lines are 
available and give consistent results, we measure extinction 
from these Balmer decrement (assuming intrinsic line ratios 



of Ha/IL3=2.87, H 7 /H/3=0.466, H<5/H/3=0.256, Osterbrock 



1989), considerably reducing the uncertainties on Ay. The 



SFR is then obtained, as in Mannucci et al. (2010), from Ha 



corrected for extinction, using the calibration in 



Kennicutt| 



(19981. Uncertainties on the SFR are computed taking into 



account the errors on both line fluxes and dust extinction. 



Finally stellar masses are taken from |Savaglio et al. ( 2009 ) . 
Table 1 lists the resulting properties of the host galaxies in 
terms of stellar mass, SFR, gas-phase metallicity, and intrin- 
sic dust extinction. 

These data are plotted in fig. [I] and compared with 
both the mass-metallicity relation (left panel) and the FMR 
(right panel) of local SDSS galaxies. We computed a lin- 
ear fit to the GRB host data taking into account the er- 
rors on metallicity, mass and SFR. The comparison with the 
mass-metallicity relation shows that, as already obtained by 
|Levesque eHST] ( |2010b| and |Han et al.| ((2010} , GRB host 
galaxies have lower metallicity than galaxies of the same 
mass both in the local universe (SDSS galaxies) and at in- 
termediate redshift ( Savaglio et al. 2005] |Jabran Zahid et"aT 



2010 ). In contrast, we also find that GRB hosts do follow the 
FMR and its extension towards low masses, without any sig- 
nificant discrepancy. In other words, when the dependence 
on SFR is properly taken into account, the metallicity prop- 
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erties of long GRB hosts do not differ substantially from 
those of the typical field population. As explained in the 
discussion, this means that the low metailicities are associ- 
ated to both low masses and high SFR, i.e, to high SSFR. 

We stress that such a good agreement is only obtained 
when the original FMR is extended using low-mass galaxies. 
The use of an extrapolation of the original 2nd-order fit 
would produce a spurious difference in metallicity, with GRB 
hosts more metal rich than field galaxies. 

In figure [2] we plot the relation between SSFR and 
metallicity for the 18 GRB host galaxies of our sample com- 
pared to the local SDSS galaxies. Here the color code shows 
different values of stellar mass. The solid lines show the re- 
lation between mass and SSFR for field galaxies, and the 
shaded area accounts for the intrinsic scatter of the observed 
relation for SDSS galaxies. GRB hosts populate the plot sim- 
ilarly to normal field galaxy population, with more (less) 
massive hosts lying in the upper (lower) bound of the ob- 
served relation. As already discussed, host metailicities are 
in line with those expected for star forming, field objects, 
apart GRB 980425. Notably, all the GRB hosts are found to 
present relatively high SSFR, with log(SSFR)^ -10. Their 
growth time, i.e. the time required by the galaxy to form 
its observed stellar mass at the present level of SFR, i.e. 
1/SSFR, is shorter than the Hubble time at the redshift of 
the GRB, for all objects in our sample. This indicates that 
GRB host are forming quite efficiently their stars similarly 
to local starbursts. 



4 DISCUSSION 

We have compared the metallicity properties of a sample 
of 18 GRB host galaxies with those of the local field pop- 
ulation. In particular, we have found that GRB hosts do 
follow the FMR recently found by |Mannucci et al.| |2010p . 
This fact implies that GRB hosts do not differ substantially 
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GRBs, and most of them reside in dusty environments (e.g., 
Perley et a l]|2009| [Kupcii Yoldas. et al.||2010[ ). It is still not 
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Figure 2. Metallicity of GRB hosts as a function of the SSFR. 
The grey-shaded area shows the relation for all SDSS galaxies, 
irrespective of mass, with the areas containing 68% and 90% of 
the galaxies. The colored lines are the median metallicity of SDSS 
galaxies with the listed values of stellar mass. Dots are GRB hosts, 
color coded with stellar mass. A broad agreement between the two 
distributions is found. 



from the typical galaxy population. The typical low, sub- 



solar metallicity found in many recent studies (e.g.,|Savaglio| 



et al. 2009 Levesque et al. 2010b and references therein) 



does not necessary mean that GRBs occur in special, low- 
metallicity galaxies, as the exception of GRB 020819 clearly 
shows, and that a direct link between low metallicity and 
GRB production exists. Indeed, this observation can be ex- 
plained as a consequence of the well-known link between 
the GRB event and the death of a very massive stars, which 
produces a relation between long GRBs and star formation 
(|Totani||1997l|Mao fc Mo||1998[|Wijers et al.|p98] |Porciani| 
|fc Madau|2001[ ). In the local universe, about 70% of all star 
formation activity o ccurs in galaxies with masses between 



lO 5 *' 5 and 10 iu / M Q (Mobasher et al. 20091, where most of 



the GRB of our sample are also found. Low stellar mass 



means low metallicities at all redshifts (e.g., |Tremonti et al 
2004||Savaglio et al.|2005||Jabran Zahid et al.|2010||Erb et aT 
2006a||Maiolino et al.|2008||Mannucci et al.|2009[ ), therefore 
low metallicities are expected for GRB hosts. Also, the FMR 
shows that galaxies with higher SFR also have lower metal- 
licities than more quiescent galaxies of the same mass. As a 
consequence, a star-formation selected galaxy sample, such 
as the present GRB host sample, is expected to fall below 
the mass-metallicity relation, but follow the FMR. This is 
what is observed in figure [I] 

Some warnings apply, which are related to the nature 
of the present GRB host galaxy sample. Our sample con- 
sists mostly of long GRBs whose position has been provided 
by the detection of their optical afterglow. It is known that 
a population of GRBs with a bright X-ray afterglow and 
lacking of optical counterpart does exist, the so-called dark 



known if this dust is spread across the host galaxy, which 
thus is likely to be metal-rich (Fynbo et al. 20091, or is 



directly associated to the GRB itself (Perley et al. 20091, 



maybe without a clear dependece on metallicity. In all cases, 
it is difficult to access with the available data whether the 
inclusion of dark GRBs could change our main conclusions. 
Indeed, it is possible that dark GRB hosts would populate 
the region of the FMR at high value of fi a . This kind of stud- 
ies will require the collection of an unbiased GRB-selected 
galaxy sample (see e.g. Malesani et al. 2009). 

Such a complete or well controlled sample is also needed 
to address the role of the several selection effects that could 
exist even within the class of GRBs with an optical bright 
afterglow. For example, it is very probable that galaxies with 
high SFR are over-represented in this sample because they 
are easier to detect. For this reason, for example, the present 
sample of GRB hosts does not allow us to study if a direct 
correspondence exist between the fraction of GRB and the 
fraction of SFR as a function of stellar mass of the hosts. 
Despite this problem, it is likely that the present sample 
represents the metallicity properties of the population of 
host of GRBs with optical bright afterglows because none 
of the main selection effects within this class of GRBs is 
directly related to metallicity. 



5 CONCLUSIONS 

Our main findings can be summarized in the following two 
points: 

(I) the average low metallicity observed in long GRB host 
galaxies is an outcome of the observed relation between stel- 
lar mass, SFR and metallicity: since the long GRB hosts are 
generally low mass and high star forming galaxies, i.e. ob- 
jects characterized by a low value of fi a , their metallicity is 
expected to be sub-solar; the metallicity observed in GRB 
hosts is exactly what is expected on the basis on their mass 
and SFR, with no apparent bias towards lower metallicities. 

(II) long GRBs not necessary explode in galaxies with a low 
metallicity (and indeed GRB 020819 is one of such cases). 
The condition for a galaxy to host a GRB seems indeed 
related to its ability to form stars in a efficient way. Our 
conclusions, based on a sample of GRB host at z<l, are 



similar to what has been found by Fynbo et al. ( 2008 1 at 
high redshift. 

Also, since GRB hosts follow the FMR, the relation can 
be used to predict the properties of those hosts for which one 
of the parameters (M*, SFR or metallicity) is not known. 

Finally, our results also suggest that larger samples 
of GRB hosts can be used to study the FMR of normal 
starburst galaxies. This is in particular exciting since GRBs 
may allow to extend current studies of the FMR both 
toward low values of fi a and toward higher redshift, in 
principle up to extremely high redshifts (at least up to 
z~8as shown by GRB 090423). Thank to they brightness, 
GRBs can be used as a background light to study the 
metal content of its parent galaxy even at very high-z as 
demonstrated by the case of GRB 050904 at z = 6.3 for 
which an estimate of the metallicity has been derived by 



Totani et al. (2006) and Kawai et al. (20061. The study 
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of the metal enrichment history at these early cosmic 
epochs is of uttermost importance to better understand 
the first stages of galaxy formation in the Universe and 
to constrain the properties of those galaxies that have re- 



ionized the intergalactic medium (see Salvaterra et al. 2010 1 
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